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Abstract 

Dial-a-Ride problems consist of a set V of n vertices in a metric space (denoting travel time 
between vertices) and a set of m objects represented as source-destination pairs {(sj, tj)}^ =1 , where 
each object requires to be moved from its source to destination vertex. In the multi-vehicle Dial-a- 
Ride problem, there are q vehicles each having capacity k and where each vehicle j € [q] has its own 
depot-vertex rj E V. A feasible schedule consists of a capacitated route for each vehicle (where 
vehicle j originates and ends at its depot rj) that together move all objects from their sources to 
destinations. The objective is to find a feasible schedule that minimizes the maximum completion 
time (i.e. makespan) of vehicles, where the completion time of vehicle j is the time when it returns 
to its depot rj at the end of its route. We study the preemptive version of multi-vehicle Dial-a-Ride, 
where an object may be left at intermediate vertices and transported by more than one vehicle, while 
being moved from source to destination. 

Our main results are an 0(log 3 n) -approximation algorithm for preemptive multi-vehicle Dial- 
a-Ride, and an improved 0(log ^-approximation for its special case when there is no capacity con- 
straint (here t < n is the number of distinct depot-vertices). There is an f^log 1 / 4-6 n) hardness 
of approximation known even for single vehicle capacitated Dial-a-Ride ifTTl . For uncapacitated 
multi-vehicle Dial-a-Ride, we show that there are instances where natural lower bounds (used in our 
algorithm) are Q(logt) factor away from the optimum. 

We also consider the special class of metrics induced by graphs excluding any fixed minor (eg. 
planar metrics). In this case, we obtain improved guarantees of 0(log 2 n) for capacitated multi- 
vehicle Dial-a-Ride, and O(l) for the uncapacitated problem. 

1 Introduction 

The multi-vehicle Dial-a-Ride problem involves routing a set of m objects from their sources to respec- 
tive destinations using a set of q vehicles starting at t distinct depot vertices in an n- vertex metric. Each 
vehicle has a capacity k which is the maximum number of objects it can carry at any time. Two versions 
arise based on whether or not the vehicle can use any vertex in the metric as a preemption (a.k.a. trans- 
shipment) point - we study the less-examined preemptive version in this paper. The objective in these 
problems is either the total completion time or the makespan (maximum completion time) over the q 
vehicles, and we study the more interesting makespan version of the problem. Thus this paper studies 
the preemptive, capacitated, minimum makespan, multi- vehicle Dial-a-Ride (mDaR) problem. 

While the multiple qualifications may make the problem appear contrived, this is exactly the problem 
that models courier or mail delivery over a day from several city depots: preemption is cheap and useful 
for packages, trucks are capacitated and the makespan reflects the daily working time limit for each 
truck. Despite its ubiquity, this problem has not been as well studied as other Dial-a-Ride versions. 
One reason from the empirical side is the difficulty in handling the possibility of preemptions in a 
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clean mathematical programming model. On the theoretical side which is the focus of this paper, the 
difficulty of using preemption in a meaningful way in an approximation algorithm persists. It is further 
compounded by the hardness of the makespan objective. 

The requirement in preemptive Dial-a-Ride, that preemptions are allowed at all vertices, may seem 
unrealistic. In practice, a subset P of the vertex-set V might represent the vertices where preemp- 
tion is permitted: the two extremes of this general problem are non-preemptive Dial-a-Ride (P = 0) 
and preemptive Dial-a-Ride (P = V). However preemptive Dial-a-Ride is more generally applica- 
ble: specifically in situations where the preemption-points P form a net of the underlying metric (i.e. 
every vertex in V has a nearby preemption-point). In particular, approximation algorithms for preemp- 
tive Dial-a-Ride imply good approximations even in this general setting, wherein the precise guarantee 
depends on how well P covers V. 

We also note that though our model allows any number of preemptions and preemptions at all ver- 
tices, our algorithms do not use this possibility to its full extent. Our algorithm for the capacitated case 
preempts each object at most once and our algorithm for the uncapacitated case only preempts objects 
at depot vertices (and at most 0(log t) times). 

The preemptive Dial-a-Ride problem has been considered earlier with a single vehicle, for which 
an O(logn) approximation [8] and an $7(log 1//4_€ n) hardness of approximation (for any constant e > 
0) ifTTl are known. Note that the completion time and makespan objectives coincide for this case. 

Moving to multiple vehicles, the total completion time objective admits a straightforward 0(log n) 
approximation along the lines of the single vehicle problem ||8l: Using the FRT tree embedding [15 ], one 
can reduce to tree-metrics at the loss of an expected 0(log n) factor, and there is a simple constant factor 
approximation for this problem on trees. The maximum completion time or makespan objective, which 
we consider in this paper turns out to be considerably harder. Due to non-linearity of the makespan 
objective, the above reduction to tree-metrics does not hold. Furthermore, the makespan objective does 
not appear easy to solve even on trees. 

Unlike in the single- vehicle case, note that an object in multi-vehicle Dial-a-Ride may be transported 
by several vehicles one after the other. Hence it is important for the vehicle routes to be coordinated so 
that the objects trace valid paths from respective sources to destinations. For example, a vehicle may 
have to wait at a vertex for other vehicles carrying common objects to arrive. The multi-vehicle Dial-a- 
Ride problem captures aspects of both machine scheduling and network design problems. 



1.1 Results, Techniques and Paper Outline 

Uncapacitated mDaR. We first consider the special case of multi- vehicle Dial-a-Ride where the vehi- 
cles have no capacity constraints (i.e. k > m). This problem is interesting in itself, and serves as a good 
starting point before we present the algorithm for the general case. The uncapacitated mDaR problem 
itself highlights differences from the single vehicle case: For example, in single vehicle Dial-a-Ride, 
preemption plays no role in the absence of capacity constraints; whereas in the multi- vehicle case, an 
optimal non-preemptive schedule may take ^(v^) longer than the optimal preemptive schedule (see 
Section[2]). Our first main result is the following: 

Theorem 1 There is an 0(log t) -approximation algorithm for uncapacitated preemptive mDaR. Addi- 
tionally, the schedule preempts objects only at depot vertices. 



The above algorithm has two main steps: the first one (in Subsection 2.1) reduces the instance at a 



constant factor loss in the performance guarantee to one in which all demands are between depots (a 



"depot-demand" instance). In the second step (Subsection 2.2 ), we use a sparse spanner on the demand 
graph to construct routes for moving objects across depots. 

We also construct instances of uncapacitated mDaR where the optimal value is Q(\ogt/ log log t) 



times all our lower bounds for this problem (Subsection 2.3 ). This suggests that stronger lower bounds 



are needed to obtain a better approximation ratio than what our approach provides. 
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We then consider the special class of metrics induced by graphs excluding some fixed minor (such 



as planar or bounded-genus graphs), and obtain the following improved guarantee in Subsection 2.4 



Theorem 2 There is an O (^-approximation algorithm for uncapacitated preemptive ITlDaR on metrics 
induced by graphs that exclude any fixed minor. 

Again the resulting schedule only preempts objects at depot-vertices. Furthermore, each object is pre- 
empted at most thrice; whereas the algorithm in Theorem [I] might preempt each object O(logt) times. 
The algorithm in Theorem[2]has the same high-level approach outlined for Theorem[T] the difference is 
in the second step, where we use a stronger notion of sparse covers in such metrics (which follows from 
the KPR decomposition theorem [22]), to construct routes for moving objects across depots. 



Capacitated mDaR. In Section [3j we study the capacitated multi-vehicle Dial-a-Ride problem, and 
obtain our second main result. Recall that there is an 0(log 1//4_e n) hardness of approximation for even 
single vehicle Dial-a-Ride ifFTl . 

Theorem 3 There is a randomized O '(log 3 n) -approximation algorithm for preemptive mDaR. Addi- 
tionally, the schedule preempts each object at most once. 

This algorithm is considerably more complex than the one for the uncapacitated special case, and is 
the main technical contribution of this paper. It has four key steps: (1) We preprocess the input so that 
demand points that are sufficiently far away from each other can be essentially decomposed into separate 
instances for the algorithm to handle independently. (2) We then solve a single- vehicle instance of the 
problem that obeys some additional bounded-delay property (Theorem [9]> that we prove; This property 
combines ideas from algorithms for light approximate shortest path trees I2TI and capacitated vehicle 
routing lfT9l . The bounded-delay property is useful in randomly partitioning the single vehicle solution 
among the q vehicles available to share this load. This random partitioning scheme is reminiscent of 
the work of Hochbaum-Maass [20], Baker [3] and Klein-Plotkin-Rao El . in trying to average out 
the effect of the cutting in the objective function. (3) The partitioned segments of the single vehicle 
tour are assigned to the available vehicles; However, to check if this assignment is feasible we solve a 
matching problem that identifies cases when this load assignment must be rebalanced. This is perhaps 
the most interesting step in the algorithm since it identifies stronger lower bounds for subproblems 
where the current load assignment is not balanced. (4) We finish up by recursing on the load rebalanced 
subproblem; An interesting feature of the recursion is that the fraction of demands that are processed 
recursively is not a fixed value (as is more common in such recursive algorithms) but is a carefully 
chosen function of the number of vehicles on which these demands have to be served. 



We prove the new bounded-delay property of s ingle -vehicle Dial-a-Ride in Subsection |3.1[ Then we 
present the algorithm for Theorempnn Subsection 3.2 and establish an 0(log 2 m log n) -approximation 



bound. Using some additional preprocessing, we show in Subsection |3.3| how to remove the dependence 
on m (number of objects) to obtain the final 0(log 3 n) approximation ratio. 

When the underlying metric is induced by graphs excluding a fixed minor, we can establish a stronger 
bounded-delay property in step (2) of the above framework. The main idea here is the construction of 
well-separated covers in such metrics, which we show can be obtained using the KPR decomposi- 



tion [22]. This leads to the following improved guarantee, proved in Subsection 3.4 

Theorem 4 There is an 0(log 2 n)-approximation algorithm for preemptive mDaR on metrics induced 
by graphs excluding any fixed minor. 

1.2 Related Work 

Dial-a-Ride problems form an interesting subclass of Vehicle Routing Problems that are well studied 
in the operations research literature. Paepe et al. lfT3l provide a classification of Dial-a-Ride problems 
using a notation similar to that for scheduling and queuing problems: preemption is one aspect in this 
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classification. Savelsberg and Sol E71 and Cordeau and Laporte lfT2l survey several variants of non- 
preemptive Dial-a-Ride problems that have been studied in the literature. Most Dial-a-Ride problems 
arising in practice involve making routing decisions for multiple vehicles. 

Dial-a-Ride problems with transshipment (the preemptive version) have been studied in Il24ll25ll26l . 
These papers consider a more general model where preemption is allowed only at a specified subset of 
vertices. Our model (and that of (H) is the special case when every vertex can serve as a preemption 
point. It is clear that preemption only reduces the cost of serving demands: Nakao and Nagamochi ||26ll 
studied the maximum decrease in the optimal cost upon introducing one preemption point. Il24ll25ll also 
model time-windows on the demands, and study heuristics and a column-generation based approach; 
they also describe applications (eg. courier service) that allow for preemptions. 

For single vehicle Dial-a-Ride, the best known approximation guarantee for the preemptive version 
is O(logn) (Charikar and Raghavachari [8]), and an Sl(log 1//4 ~ e n) hardness of approximation (for any 
constant e > 0) is shown in G0rtz ifTTl . The non-preemptive version appears much harder and the 
best known approximation ratio is Tam.{^/k log n, ^fn log 2 n} (Charikar and Raghavachari [ 8 ], Gupta et 
al. ifTHl ): however to the best of our knowledge, APX-hardness is the best lower bound. There are known 
instances of single vehicle Dial-a-Ride where the ratio between optimal non-preemptive and preemptive 
tours is Q,{y/n) in general metrics JH, and ^(n 1 / 8 ) in the Euclidean plane ifTBl . A 1.8-approximation is 
known for the k = 1 special case of single vehicle Dial-a-Ride (a.k.a. stacker-crane problem) |fl6l . 

Charikar et al. Q studied the related k-delivery TSP problem, which involves transporting a number 
of identical objects from supply to demand vertices, using a single capacity k vehicle. The key difference 
from Dial-a-Ride is that an object can be moved from any supply vertex to any demand vertex. [7] gave 
an approximation algorithm for fc-delivery TSP that outputs a non-preemptive tour of length at most five 
times an optimal preemptive tour. They also showed that for any fc-delivery TSP instance, the optimal 
non-preemptive tour is at most four times the optimal preemptive tour. 

The truck and trailer routing problem J5j [28j is another problem where preemption plays a crucial 
role. Here a number of capacitated trucks and trailers are used to deliver all objects. Some customers 
are only accessible without the trailer. The trailers can be parked at any point accessible with a trailer 
and it is possible to shift demand loads between the truck and the trailer at the parking places. The 
papers [51 |28l present heuristics for this problem. 

Single vehicle preemptive Dial-a-Ride is closely related to the uniform buy-at-bulk problem (with 
cost function [f ] where k is the vehicle capacity). Such a connection was formally used in G0rtz IfTTl to 
establish the hardness of approximation for the single vehicle problem. Approximation algorithms for 
several buy-at-bulk variants have been studied recently, eg. non-uniform buy-at-bulk |6l|9], buy-at-bulk 
with node-costs iflOll and buy-at-bulk with protection Q ; poly-logarithmic approximation guarantees are 
known for all these problems. However the techniques required to solve the multi-vehicle Dial-a-Ride 
problem appear quite different from these buy-at-bulk results. 

The uncapacitated imDaR problem generalizes the nurse-station-location problem that was studied 
in Even et al. lfl4l (where a 4-approximation algorithm was given). In fact we also use this algorithm as a 
subroutine for uncapacitated mDaR. Nurse-station-location is the special case of uncapacitated mDaR 
when each source-destination pair coincides on a single vertex. In this paper, we handle not only the 
case with arbitrary pairs (uncapacitated mDaR), but also the more general problem (capacitated mDaR) 
with finite capacity restriction. 

1.3 Problem Definition and Preliminaries 

We represent a finite metric as (V, d) where V is the set of vertices and d is a symmetric distance function 
satisfying the triangle inequality. For subsets A, B C V we denote by d(A, B) the minimum distance 
between a vertex in A and another in B, so d(A, B) = mm{d(u, v) \ u 6 A, v G B}. For a subset 
E C (^) of edges, d(E) := J2eeE de denotes the total length of edges in E. 

The multi-vehicle Dial-a-Ride problem (mDaR) consists of an n-vertex metric (V, d), m objects 
specified as source-destination pairs {s^ti}™!, q vehicles having respective depot-vertices {Vj}| =1 , 
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and a common vehicle capacity k. A feasible schedule is a set of q routes, one for each vehicle (where 
the route for vehicle j G [q] starts and ends at rj), such that no vehicle carries more than k objects 
at any time and each object is moved from its source to destination. The completion time Cj of any 
vehicle j G [q] is the time when vehicle j returns to its depot rj at the end of its route (the schedule is 
assumed to start at time 0). The objective in mDaR is to minimize the makespan, i.e., minmax Jg [ 9 ] Cj. 
We denote by S := {si \ i G [to]} the set of sources, T := {t{ \ i G [to]} the set of destinations, 
R := {rj | j G [q]} the set of distinct depot- vertices, and t := \R\ the number of distinct depots. In this 
paper, we only consider the preemptive version, where objects may be left at intermediate vertices while 
being moved from source to destination. 

Single vehicle Dial-a-Ride. The following are lower bounds for the single vehicle problem: the 
minimum length TSP tour on the depot and all source/destination vertices (Steiner lower bound), and 
i=1 k - (flow lower bound). Charikar and Raghavachari [8] gave an 0(log n)-approximation algo- 
rithm for this problem based on the above lower bounds. A feasible solution to preemptive Dial-a-Ride 
is said to be 1 -preemptive if every object is preempted at most once while being moved from its source 
to destination. Gupta et al. [ 18 ] showed that the single vehicle preemptive Dial-a-Ride problem always 
has a 1-preemptive tour of length 0(log 2 n) times the Steiner and flow lower-bounds. 

y^ m d(s t ) 

Lower bounds for mDaR. The quantity i=1 fc " is a lower bound similar to the flow bound for 
single vehicle Dial-a-Ride. Analogous to the Steiner lower bound above, is the optimal value of an in- 
duced nurse-station-location instance. In the nurse-station-location problem [14), we are given a metric 
(V,d), a set T of terminals and a multi-set {r 3 }| =1 of depot- vertices; the goal is to find a collection 
{Fj} q - =l of trees that collectively contain all terminals T such that each tree Fj is rooted at vertex rj and 
max j=i d(Fj) is minimized. Even et al. |[T4l gave a 4-approximation algorithm for this problem. The 
optimal value of the nurse-station-location instance with depots {rj} q - =l (depots of vehicles in mDaR) 
and terminals T = S U T is a lower bound for mDaR. The following are some lower bounds implied 
by nurse-station-location: (a) 1/q times the minimum length forest that connects every vertex in S U T 
to some depot vertex, (b) max ig [ m ] d(R, Si), and (c) maxj g r m i d(R, ti). Finally, it is easy to see that 
maxj g r m i d(si,ti) is also a lower bound for mDaR. 

We note that our approximation bounds for uncapacitated mDaR are relative to the above lower 
bounds. However the algorithm for capacitated mDaR relies additionally on stronger lower bounds 
derived from suitable subproblems. 



2 Uncapacitated Multi- Vehicle Dial-a-Ride 

In this section we study the uncapacitated special case of mDaR, where vehicles have no capacity con- 
straints (i.e. capacity k > to). We give an algorithm that achieves an O(logt) approximation ratio for 
this problem (recall t < n is the number of distinct depots). Unlike in the single vehicle case, preemptive 
and non-preemptive versions of mDaR are very different even without capacity constraints. 

Preemption gap in Uncapacitated mDaR. Consider an instance of uncapacitated mDaR where the 
metric is induced by an unweighted star with q leaves (where q is number of vehicles), all q vehicles 
have the center vertex as depot, and there is a demand between every pair of leaf-vertices. A preemptive 
schedule having makespan 4 is as follows: each vehicle j G [q] visits leaf j and brings all demands 
with source j to the root, then each vehicle j visits its corresponding leaf again, this time delivering all 
demands with destination j. On the other hand, in any non-preemptive schedule, one of the q vehicles 
completely serves at least q — 1 demands (since there are q(q — 1) demands in all). The minimum 
length of any tour containing the end points of q demands is Q (y/q), which is also a lower bound on the 
optimal non-preemptive makespan. Thus there is an ^(y^) factor gap between optimal preemptive and 
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non-preemptive tours. This is in contrast to the uncapacitated single vehicle case, where it is easy to see 
that the optimal preemptive and non-preemptive tours coincide. 

The algorithm for uncapacitated imDaR proceeds in two stages. Given any instance, it is first reduced 
(at the loss of a constant factor) to a depot-demand instance, where all demands are between depot 



vertices (Subsection 2.1 1. This reduction uses the nurse-station-location algorithm from Even et al. 111411 . 



Then the depot-demand instance is solved using an 0(log ^-approximation algorithm (Subsection 2.2 1; 
this is the key step in the algorithm, and is based on a sparse-spanner on the demand graph. 

2.1 Reduction to Depot-demand Instances 

We define depot-demand instances as those instances of uncapacitated mDaR where all demands are 
between depot vertices. Given any instance X of uncapacitated mDaR, the algorithm UncapMuIti (Fig- 
ure[T]) reduces X to a depot-demand instance. We now argue that the reduction in UncapMuIti only loses 
a constant approximation factor. Let B denote the optimal makespan of instance X. Since the optimal 
value of the nurse- station-location instance solved in the first step of UncapMuIti is a lower bound for 
X, we have max' =1 d(Fj) < AB. 



Input: instance I of uncapacitated mDaR. 

• Solve the nurse-station-location instance with depots {rj}j =1 and all sources/destinations 
S U T as terminals, using the 4-approximation algorithm fPfl . Let {-F)}?=i be the resulting 
trees covering S UT such that each tree Fj is rooted at depot rj. 

• Define a depot-demand instance J of uncapacitated mDaR on the same metric and set of 
vehicles, where the demands are {(rj, rj) | s$ € Fj & ti € Fi, 1 < i < m}. For any object 
i E [m] let the source depot be the depot rj for which Sj <E Fj and the destination depot be 
the depot r; for which f, e F|. 

• Output the following schedule for I: 

1 . Each vehicle j e [q] traverses tree Fj by an Euler tour, picks up all objects from sources 
in Fj and brings them to their source-depot Tj . 

2. Vehicles implement a schedule for depot-demand instance J , and all objects are moved 



from their source-depot to destination-depot (using the algorithm in Section 2.2 1. 

3. Each vehicle j G [q] traverses tree Fj by an Euler tour, picks up all objects having 
destination-depot Tj and brings them to their destinations in Fj. 

Figure 1: Algorithm UncapMuIti for uncapacitated mDaR. 



Claim 5 The optimal makespan for the depot-demand instance J is at most 17 B. 

Proof: Consider a feasible schedule for J involving three rounds: (1) each vehicle traverses (by means 
of an Euler tour) its corresponding tree in {Fj} q - =1 and moves each object i from its source-depot (the 
source in instance J) to Si (source in original instance X); (2) each vehicle follows the optimal schedule 
for X and moves each object % from S{ to ti (destination in X); (3) each vehicle traverses its corresponding 
tree in {Fj} q - =l and moves each object i from U to its destination-depot (the destination in J). Clearly 
this is a feasible schedule for J . From the observation on the nurse-station-location instance, the time 
taken in each of the first and third rounds is at most 8B. Furthermore, the time taken in the second round 
is the optimal makespan of X which is B. This proves the claim. ■ 

Assuming a feasible schedule for J, it is clear that the schedule returned by UncapMuIti is feasible 
for the original instance X. The first and third rounds in X's schedule require at most 8B time each. 
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Thus an approximation ratio a for depot-demand instances implies an approximation ratio of 17a + 8 
for general instances. In the next subsection, we show an O (log t) -approximation algorithm for depot- 
demand instances (here t is the number of depots), which implies Theorem[T] 

2.2 Algorithm for Depot-demand Instances 

Let J be any depot-demand instance: note that the instance defined in the second step of UncapMuIti is 
of this form. It suffices to restrict the algorithm to the induced metric (R, d) on only depot vertices, and 
use only one vehicle at each depot in R. Consider an undirected graph H consisting of vertex set R and 
edges corresponding to demands: there is an edge between vertices r and s iff there is an object going 
from either r to s or s to r. Note that the metric length of any edge in H is at most the optimal makespan 
B of instance J. In the schedule produced by our algorithm, vehicles will only use edges of H. Thus in 
order to obtain an O(logi) approximation, it suffices to show that each vehicle only traverses O(logt) 
edges. Based on this, we further reduce J to the following instance T~L of uncapacitated imDaR: the 
underlying metric is shortest paths in graph H (on vertices R), with one vehicle at each i?-vertex, and 
for every edge (u, v) G H there is a demand from u to v and one from v to u. Clearly any schedule for 
7-L having makespan /3 implies one for J of makespan (3 ■ B. The next lemma implies an 0(log \R\) 
approximation for depot-demand instances. 

Lemma 6 There exists a poly-time computable schedule for 7i with makespan 0(log t), where t = \R\. 

Proof: Let a = [~lg t] + 1. We first construct a sparse spanner A of H as follows: consider edges of 
H in an arbitrary order, and add an edge (u, v) G H to A iff the shortest path between u and v using 
current edges of A is more than 2a. It is clear from this construction that the girth of A (length of its 
shortest cycle) is at least 2a, and that for every edge (u, v) G H, the shortest path between u and v in A 
is at most 2a. 

We now assign each edge of A to one of its end-points such that each vertex is assigned at most 
two edges. Repeatedly pick any vertex v of degree at most two in A, assign its adjacent edges to v, 
and remove these edges and v from A. We claim that at the end of this procedure (when no vertex has 
degree at most 2), all edges of A would have been removed (i.e. assigned to some vertex). Suppose for a 
contradiction that this is not the case. Let A ^ be the remaining graph; note that A C A, so girth of A 
is at least 2a. Every vertex in A has degree at least 3, and there is at least one such vertex w. Consider 
performing a breadth-first search in A from w. Since the girth of A is at least 2a, the first a levels of 
the breadth-first search is a tree. Furthermore every vertex has degree at least 3, so each vertex in the 
first a — 1 levels has at least 2 children. This implies that A has at least 1 + 2 a ~ 1 > t vertices, which 
is a contradiction! For each vertex v G R, let A v denote the edges of A assigned to v by the above 
procedure; we argued that U v& rA v = A, and \A V \ < 2 for all v G R. 

The schedule for 7-L involves 2a rounds as follows. In each round, every vehicle v G R traverses the 
edges in A v (in both directions) and returns to v. Since \A V \ < 2 for all vertices v, each round takes 4 
units of time; so the makespan of this schedule is 8a = 0(log t). The route followed by each object in 
this schedule is the shortest path from its source to destination in spanner A; note that the length of any 
such path is at most 2a. To see that this is indeed feasible, observe that every edge of A is traversed by 
some vehicle in each round. Hence in each round, every object traverses one edge along its shortest path 
(unless it is already at its destination). Thus after 2a rounds, all objects are at their destinations. ■ 

Note that the final algorithm for uncapacitated mDaR preempts each object at most O(logt) times, 
and only at depot-vertices. This completes the proof of Theorem[T] 

2.3 Tight Example for Uncapacitated mDaR Lower Bounds 

We note that known lower bounds for uncapacitated mDaR are insufficient to obtain a sub-logarithmic 
approximation guarantee. The lower bounds we used in our algorithm are the following: max ie [ m ] d(si,ti), 
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and the optimal value of a nurse-station- location instance with depots {rj}J =1 and terminals SUT. We 
are not aware of any lower bounds stronger than these two bounds. 

We show an instance Q of uncapacitated mDaR where the optimal makespan is a factor 
larger than both the above lower bounds. In fact, the instance we construct is a depot-demand instance 
that has the same special structure as instance % in Lemma [6] I.e. the demand graph is same as the 
graph inducing distances. Take G = (V, E) to be a i-vertex regular graph of degree ~ log t and girth 
g ~ logt/ loglogt (there exist such graphs, eg. Lazebnik et al. (23l). Instance Q is defined on a metric 
on vertices V with distances being shortest paths in graph G. For every edge (u, v ) G E of graph G, 
there is an object with source u and destination v (the direction is arbitrary). There is one vehicle located 
at every vertex of V; so number of vehicles q = t. 

Observe that both our lower bounds are O(l): the optimal value of the nurse-station- location in- 
stance is 0, and maximum source-destination distance is 1. However as we show below, the optimal 
makespan for this instance is at least g — 1 = f2(log tj loglog t). Suppose (for contradiction) that there 
is a feasible schedule for Q with makespan M < g — 1. A demand (u, v ) G E is said to be completely 
served by a vehicle j iff the route of vehicle j visits both vertices u and v. The number of distinct 
vertices visited by any single vehicle is at most M < g: so the number of demands that are completely 
served by a single vehicle is at most M — 1 (otherwise these demand edges would induce a cycle smaller 
than the girth g). Hence the number of demands that are completely served by some vehicle is at most 
t • g < \E\. Let (u, v) G E be a demand that is not completely served by any vehicle, i.e. there is no 
vehicle that visits both u and v. Since we have a feasible schedule of makespan M, the path ir followed 
by demand (u, v) from u to v (or vice versa) in the schedule has length at most M. The path it can not 
be the direct edge (u, v) since demand (u, v) is not completely served by any vehicle. So path it together 
with edge (u, v) is a cycle of length at most M + 1 < g in graph G, contradicting girth of G. 

2.4 Improved Algorithm for Metrics Excluding a Fixed Minor 

We now prove Theorem [2] that gives a constant approximation algorithm for uncapacitated mDaR on 
metrics induced by i^-minor free graphs (for any fixed r). This improvement comes from using the 
existence of good 'sparse covers' in such metrics, as opposed to the spanner based construction in 
Lemma[6] This guarantee is again relative to the above mentioned lower bounds. 

Consider an instance of uncapacitated mDaR on metric (V, d) that is induced by an edge- weighted 
graph G = (V, E) containing no iv~ r -minor (for some fixed r > 1). We start with some definitions @. 
A cluster is any subset of vertices. For any 7 > and vertex v G V, N(v, 7) := {u G V \ d(u, v) < 7} 
denotes the set of vertices within distance 7 from v. As observed in Busch et al. [4J and Abraham et 
al. HI, the partitioning scheme of Klein et al. [22] implies the following result. 

Theorem 7 (|22|) Given K r -minor free graph G = (V, E, w) and value 7 > 0, there is an algorithm 
that computes a set Z = {C\, ■ ■ ■ , C{\ of clusters satisfying: 

1. The diameter of each cluster is at most 0(r 2 ) • 7, i.e. max u v ^d d(u, v) < 0(r 2 ) • 7/or all i G [I]. 

2. For every v G V, there is some cluster d G Z such that N(v, 7) C Cj. 

3. For every v G V, the number of clusters in Z that contain v is at most 0{2 r ). 
The set of clusters Z found above is called a sparse cover of G. 

Proof: [Theorem |2) The reduction in Section [2] implies that it suffices to consider depot-demand in- 
stances: An O(l) approximation for such instances implies an O(l) approximation for general instances. 
Let J be any depot-demand instance on metric (V, d) induced by K T -minor free graph G = (V,E), with 
a set R C V of depot-vertices (each containing a vehicle), and where all demands {si,ti} 7 ^ =l are between 
vertices of R. The algorithm is described in Figure [2] 
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Input: Depot-demand instance J on metric G = (V, E, w), depot-vertices R C V, demands {s.;, 

1. Let 7 = max ie [ m ] d(si,ti) be the maximum source-destination distance. 

2. Compute a sparse cover Z = {Cj} l j =1 given in Theoremj7]for parameter 7. 

3. For each cluster Cj € Z, choose an arbitrary vertex cj G Cj as its center. 

4. For each demand i € [to], let 7r(i) € [Z] be such that Sj, ^ £ C n (i). 

5. Output the following schedule for J: 

(a) Each vehicle r £ R visits the centers of all clusters containing r, and returns to r. Vehicle r 
carries all the objects {i e [to] \ s t — r} having source r, and drops each object i at center 

Cjr(i)- 

(b) Each vehicle r E R again visits the centers of all clusters containing r. Vehicle r brings all 
the objects {i G [to] | ti = r} having destination r: each object i is picked up from center 

c 7r ( l )- 

Output: An 0(l)-approximate minimum makespan schedule for J . 



Figure 2: Algorithm for uncapacitated mDaR on i^-minor free graphs. 



Note that 7 is a lower bound on the optimal makespan of J . We claim that the makespan of the 
above schedule is at most 0(r 2 2 r ) • 7. Observe that each depot is contained in at most 0(2 r ) clusters, 
and the distance from any depot to the center of any cluster containing it is at most 0(r 2 ) • 7. Hence the 
time taken by each vehicle in either of the two rounds (Steps (|5a|-(5b 1) is at most 0(r 2 2 r ) ■ 7. Since r 
is a fixed constant, the final makespan is O(l) -7. 

We now argue the feasibility of the above schedule. Step Q is well-defined: for all i G [m], we 
have Si,ti G N(si, 7) and there is some j G [I] with N(si, 7) C Cj (i.e. we can set = j). It is now 
easy to see that each object % G [m] traces the following route in the final schedule: S{ ~~> c„u) tj. ■ 



we 



Combining this algorithm for depot-demand instances with the reduction in Subsection 2.1 
obtain an 0(r 2 2 r ) -approximation algorithm for uncapacitated mDaR on K r -minor free metrics. When 
r is a fixed constant this is a constant approximation, and we obtain Theorem [2] Note that each object 
in the resulting schedule is preempted at most thrice, and only at depot- vertices. 



3 Capacitated Multi-Vehicle Dial-a-Ride 



In this section we prove our main result: an 0(log 3 n)-approximation algorithm for the mDaR problem 
(with capacities). We first obtain a new structure theorem on single-vehicle Dial-a-Ride tours (Subsec- 
tion |3.1| > that preempts each object at most once, and where the total time spent by objects in the vehicle 
is bounded. Obtai ning such a single vehicle tour is crucial in our algorithm for capacitated mDaR, which 
app ears in Section 3.2 There we prove a weaker approximation bound of 0(log 2 m logn); in Subsec- 
tion 3.3 we show how to remove the dependence on m to obtain an 0(log 3 n) -approximation algorithm 
even for "weighted mDaR" (i.e. Theorem [3 1. In Subsection |3.4| we establish a stronger structure the- 
orem (analogous to the one in Subsection |3.1[ ) for metrics excluding any fixed minor; this immediately 
leads to an improved 0(log 2 n)-approximation ratio for mDaR on such metrics (i.e. Theorem]^. 



3.1 Capacitated Vehicle Routing with Bounded Delay 

Before we present the structural result on Dial-a-Ride tours, we consider the classic capacitated vehicle 
routing problem |fT9l with an additional constraint on object 'delays'. The capacitated vehicle routing 
problem (CVRP) is a special case of single vehicle Dial-a-Ride when all objects have the same source (or 
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equivalently, same destination). Formally, we are given a metric (V, d), specified depot- vertex r £ V, 
and m objects all having source r and respective destinations {£i}j 6 r m ]. The goal is to compute a 
minimum length non-preemptive tour of a capacity k vehicle originating at r that moves all objects from 
r to their destinations. In CVRP with bounded delay, we are additionally given a delay parameter (3 > 1, 
and the goal is to find a minimum length capacitated non-preemptive tour serving all objects such that 
the time spent by each object i G [m] in the vehicle is at most j3 • d(r, U). Well-known lower bounds 
for the CVRP JT9ll are as follows: minimum length TSP tour on {r} U {t{ \ i G [m]} (Steiner lower 
bound), and | YliLi d{r, ti) (flow lower bound). These lower bounds also hold for the (less constrained) 
preemptive version of CVRP. 

Theorem 8 There is a (2.5 + g^j) -approximation algorithm for CVRP with bounded delay, where 
(3 > 1 is the delay parameter. This guarantee is relative to the Steiner and flow lower bounds. 

Proof: Our algorithm is basically a combination of the algorithms for light approximate shortest path 
trees (2TJ, and capacitated vehicle routing |fl9l . Let LB denote the maximum of the Steiner and flow 
lower bounds. The minimum TSP tour on the destinations plus r is the Steiner lower bound. The first 
step is to compute an approximately minimum TSP tour C: Christofides' algorithm [11 J gives a 1.5- 
approximation, so d(C) < 1.5 • LB. Number the vertices in the order in which they appear in C, starting 
with r being 0. Using the procedure in Khuller et al. J2TJ, we obtain a set of edges {(0, vi), • • • , (0, vt)} 
with < v\ < V2 < • • • < vt < \V\ having the following properties (below vo = 0). 

1. For 1 < p < t, for any vertex u with v p -\ < u < v p , the length of edge (0, v p -\) plus the path 
along C from u p _i to u is at most f3 ■ d(0, u). 

2- J2Ud(0,v p )<^-d(C). 

For each 1 < p < t, define tour C p which starts at r, goes to traverses C until v p , then returns to 
r. Assign vertices {v p -i, • • • , v p — 1} (and all demands contained in them) to C p . Also define tour Ct+\ 
which starts at r, goes to vt, and traverses C until r; and assign all remaining demands to Ct+\- It is clear 
that C p (for 1 < p < i+ 1) visits each vertex assigned to it within j3 times the shortest path from r (using 
property 1 above). Also, the total length Yl P t\ d(C p ) = d{C) + 2 £* =1 d(0, v p ) < (1 + ^ I )d(C), by 
property 2. 

For each C p , we service the set D p of demands assigned to it separately. Index the demands in D p 
in the order in which they appear on C p (breaking ties arbitrarily). Consider a capacitated tour which 
serves these demands in groups of at most k each, and returns to r after serving each group. The groups 
are defined as follows: starting at index 1, each group contains the next k contiguous demands (until all 
of D p is assigned to groups). By rotating the indexing of demands, there are k different groupings of D p 
that can be obtained: each of which corresponds to a capacitated tour. As argued in [ 19 ] (and is easy to 
see), the average length of these k tours is at most d{C p ) + 2 ^2 zeDp d ^ z ^ . So the minimum length tour 

7 P among these satisfies d(^ p ) < d(C p ) + 2 Y^zeD 

The final solution 7 is the concatenation of tours 71 , • • • , 7t+i - Note that the time spent in the vehicle 
by any demand % is at most f3 ■ d(r, ti). The length of tour 7 is at most X^p=i d(C p ) + 2 ^2 zeD » 

where D is the set of all demands. Note that 2 J2 zeD is the flow lower bound for this CVRP 
instance, so it is at most LB. Hence we obtain the following. 

d{i) < (1 + J^iWC) + LB < (1 + ^ rI )^LB + LB = (2.5 + t^t) 16 

Clearly, solution 7 satisfies the desired conditions in the theorem. ■ 

We now consider the single vehicle preemptive Dial-a-Ride problem given by a metric (V, d), a set 
D of demand-pairs, and a vehicle of capacity k. Given Theorem [8} the following structural result is a 
simple extension of Theorem 16 from Gupta et al. [ 18 ]. We present the proof for completeness. 
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Theorem 9 There is a randomized poly-time computable 1 -preemptive tour r servicing D that satisfies 
the following conditions (where LB is the average of the Steiner and flow lower bounds): 

1. Total length: d(r) < 0(log 2 n) • LB. 

2. Bounded delay: J2ieD ^ — 0(log n) J2ieD d(si,U) where Tj is the total time spent by object 
i € D in the vehicle under the schedule given by t. 

Proof: We follow the algorithm of Gupta et al. [18] that obtains a single vehicle 1 -preemptive tour 
within 0(log 2 ra) factor of the Steiner and flow lower bounds. In addition, we will also ensure the 
bounded-delay property. Using the results on probabilistic tree embedding |fl5h we may assume that 
the given metric is a hierarchically well-separated tree T with distance function n. This only increases 
the expected value of the lower bound by a factor of O(logra); i.e. £J[LB] = O(logn) LB where LB 
(resp. LB) equals the average of the Steiner and flow lower bounds on metric k (resp. d). Furthermore, 
as discussed in [18], we may assume that the metric k is induced by a tree T on the original vertex set 
V having I = O(logra) levels. 

We now partition the demands in T into I sets with D p (for p = 1, ■ • • , I) consisting of all demands 
% € D having their nearest common ancestor (nca) in level p (i.e. the nca of Sj and ti is a vertex in level 
p). We service each D p separately using a tour of length O(LB). Finally we concatenate the tours for 
each level p, to obtain the theorem. 

Servicing D p . For each vertex v at level p in T, let L v denote the demands in D p that have v as their nca. 
Let LB(<y) denote the average of the Steiner and flow lower-bounds (in metric k) for the single vehicle 
problem with depot v and demands L v . Clearly the flow lower bounds are disjoint for different vertices 
v. Also, the subtrees under any two different level p vertices are disjoint, so the Steiner lower bounds 
are disjoint for different v. Thus J2 V LB(u) < LB. We now show how each L v is served separately. 

Servicing L v . Consider the following two instances of the CVRP problem in metric k: T src with all 
sources in L v and common destination v, and Td es t with common source v and all destinations in L v . 
Observe that the Steiner and flow lower bounds for each of I src and Idest are all at most 2 LB V . 

Consider instance T src . Setting delay parameter f3 = 2 in Theorem [8} we obtain a non-preemptive 
tour a v that moves all objects in L v from their sources to vertex v, such that (1) the length of a v is at 
most 11 LB V , and (2) the time spent by each object i G L v in the vehicle is at most 2 n(si, v). Similarly 
for instance Idest, we obtain a non-preemptive tour t v that moves all objects in L v from vertex v to 
their destinations, such that (1) the length of r v is at most 11 LBt,, and (2) the time spent by each object 
i G L v in the vehicle is at most 2 k(v, ti). Concatenating these two tours, we obtain that a v ■ t v is a 
1-preemptive tour servicing L v , of length at most 22 • LB(v), where each object i £ L v spends at most 
2(k(sj, v) + k(v, ti)) = 2 • K(si,ti) time in the vehicle. 

We now obtain a depth-frrst-search traversal (DFS) on T (restricted to the end-points of demands D p ) to 
visit all vertices v in level p that have some demand in their subtree (i.e. L v ^ 0), and use the algorithm 
described above for servicing demands L v when v is visited in the DFS. This is the tour servicing D p . 
Note that the length of the DFS is at most the Steiner lower bound on k, so it is at most 2 LB. Thus the 
tour servicing D p has length at most 2 LB + 22 J2 V LB(u), where v ranges over all vertices in level p of 
T. Recall that Yl v LB(u) < LB, so the tour servicing D p has length at most 24 • LB. Additionally, the 
time spent by each object i e D p in the vehicle is at most 2 K(sj, t{). 

Finally concatenating the tours for each level p = 1, • • • ,1, we obtain a 1-preemptive tour on T of 
length O(logra) • LB. Additionally, the time spent by each object i in the vehicle is at most 2 • n(si, ti). 
This translates to a 1-preemptive tour on the original metric having expected length 0(log 2 n) ■ LB, and 
with E\%2 ieD Ti] < J2ieD 2 ' ^ U)] < 0(log n) J2ieD d ( s ii U), where T, is the time spent by 
object i in the vehicle. Using Markov inequality and a union bound, we can ensure that the resulting 
1-preemptive tour on metric d satisfies both conditions 1-2 with constant probability. ■ 
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3.2 Algorithm for Capacitated mDaR 



We are now ready to present our algorithm for capacitated multi-vehicle Dial-a-Ride. The algorithm for 
mDaR relies on a partial coverage algorithm Partial that given subsets Q C [q] of vehicles and D C [m] 
of demands, outputs a schedule for Q of near-optimal makespan that covers some fraction of demands 
in D. The main steps in Partial are as follows. 

1. Obtain a single-vehicle tour satisfying 1-preemptive and bounded-delay properties (Theorem[9]>. 

2. Randomly partition the single vehicle tour into \Q\ equally spaced pieces. 

3. Solve a matching problem to assign some of these pieces to vehicles of Q, so as to satisfy a subset 
of demands in D. 



4. A suitable fraction of unsatisfied demands in D are covered recursively by unused vehicles of Q. 

The algorithm first guesses the optimal makespan B of the given instance of mDaR (it suffices to 
know B within a constant factor, which is required for a polynomial-time algorithm). Let parameter 
1 - j+jL^. For any subset PC [q], we abuse notation and use P to denote both the set of vehicles 
P and the multi-set of depots corresponding to vehicles P. 

We give an algorithm Partial that takes as input a tuple (Q, D, B) where Q C [q] is a subset of 
vehicles, D C [m] a subset of demands and B € M. + , with the promise that vehicles Q (originating 
at their respective depots) suffice to completely serve the demands D at a makespan of B. Given such 
a promise, Partial {Q,D,B} returns a schedule of makespan O(logrelogm) • B that serves a good 
fraction of D. Algorithm Partial(Q, D, B) is given in Figure [3] We set parameter p = 6 (log ra log to), 
the precise constant in the ©-notation comes from the analysis. 



Lemma 10 If there exists a schedule of vehicles Q covering all demands D, having makespan at most 
B, then Partial invoked on (Q, D, B) returns a schedule of vehicles Q of makespan at most (16+ 16/)) B 
that covers at least an a lgz fraction of D, where z := min{|Q|, 2m} < 2m. 



The final algorithm invokes Partial iteratively until all demands are covered: each time with the entire set 
[q] of vehicles, all uncovered demands, and bound B. If D C [to] is the set of uncovered demands at any 
iteration, Lemma [To| implies that Partial ([q], D, B) returns a schedule of makespan 0(logmlog n) • B 
that serves at least?] JD| demands. Hence a standard set-cover analysis implies that all demands will be 
covered in 0(log m) rounds, resulting in a makespan of 0(log 2 to log n) ■ B. 



Proof of Lemma 10 We proceed by induction on the number \Q\ of vehicles. The base case \Q\ = 1 is 
easy: the tour r in Step ([2]) has length 0(log 2 n) ■ B < pB, and satisfies all demands (i.e. fraction 1). In 
the following, we prove the inductive step, when \Q\ > 2. 



Preprocessing. Suppose Step ([I]) applies. Note that d{V\, V2) > B and hence there is no demand with 
source in one of {Vi, V2} and destination in the other. So demands D\ and D2 partition D. Furthermore 
in the optimal schedule, vehicles Qj (any j = 1,2) only visit vertices in Vj (otherwise the makespan 
would be greater than B). Thus the two recursive calls to Partial satisfy the assumption: there is some 
schedule of vehicles Qj serving Dj having makespan B. Inductively, the schedule returned by Partial 
for each j = 1, 2 has makespan at most (16 + 16p) ■ B and covers at least a lgc • \Dj\ demands from 
Dj, where c < min{|<5| — 1, 2m} < z. The schedules returned by the two recursive calls to Partial can 
clearly be run in parallel and this covers at least a lg ^(|-Di| + I-D2I) demands, i.e. an a lgz fraction of D. 
So we have the desired performance in this case. 
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Input: Vehicles Q C [q], demands D C [m] , bound B > such that Q can serve all demands in D 
at makespan B. 

Preprocessing 

1. If the minimum spanning tree (MST) on vertices Q contains an edge of length greater than 
3B, there is a non-trivial partition {Qi, Q2} of Q with d(Q±, Q2) > SB. For j G {1, 2}, let 
Vj ■ = {v G V I d(Qj,v) < B} and be all demands of Z) induced on Vj. Run in parallel the 
schedules from Partial(Qi, D\, B) and Partial {Q25 D2, B). Assume this is not the case in the 
following. 

Random partitioning 

2. Obtain single-vehicle 1-preemptive tour r using capacity k and serving demands D, by apply- 
ing Theorem [9] 

3. Choose a uniformly random offset rj G [0, pB] and cut edges of tour r at distances {p/)-B + 77 | 
p = 1, 2, • • • } along the tour to obtain a set V of pieces of r. 

4. C" is the set of objects i G D such that 2 is carried by the vehicle in r over some edge that is 
cut in Step ((3]); and C := D\ C". Ignore C" objects in the rest of the algorithm. 

Load rebalancing 

5. Construct bipartite graph H with vertex sets V and Q and an edge between piece P G V and 
depot / G Q iff d(f, P) < 2B. For any subset 4Cp, T(yl) C Q denotes the neighborhood 

I cl 

of A in graph H. Let S C be any maximal set that satisfies < 

6. Compute a 2-matching it :V \ S —> Q\ T(S), i.e. function s.t. (P, 7r(P)) is an edge in if for 
all P G P \ S, and the number of pieces mapping to any / G Q\ T(S) is |vr _1 (/)| < 2. 

Recursion 

7. Define G\ := {i G C" | either Si G 5 or U G 5}; and C 2 := C'\Ci. 

8. Run in parallel the recursive schedule Partial(r(<S), C\,B) for Ci and the following for C2: 

(a) Each vehicle / G Q \ r(<S) traverses the pieces 7r _1 (/), moving all ^-objects in them 
from their source to preemption- vertex, and returns to its depot. 

(b) Each vehicle / G Q \ T(S) again traverses the pieces 7r _1 (/), this time moving all C2- 
objects in them from their preemption- vertex to destination, and returns to its depot. 

Output: A schedule of Q of makespan (16 + 16p) • B that serves an ct lgmm ^^'' 2m ^ fraction of D. 



Figure 3: Algorithm Partial(Q, D, B) for capacitated mDaR. 



Random partitioning. The harder part of the algorithm is when Step ([T]) does not apply: so the MST 
length on Q is at most 3|Q| ■ B. Note that when the depots Q are contracted to a single vertex, the MST 
on the end-points of D plus the contracted depot- vertex has length at most \Q\-B (the optimal makespan 
schedule induces such a tree). Thus the MST on the depots Q along with end-points of D has length at 
most 4\Q\ • B. Based on the assumption in Lemma 10 and the flow lower bound for mDaR, we have 
Yli^D d(si,U) < k\Q\ ■ B. It follows that for the single vehicle Dial-a-Ride instance solved in Step Q, 
the Steiner and flow lower-bounds (denoted LB in Theorem[9]) are 0(1) • \Q\B. Theorem[9]now implies 
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that (with high probability) r is a 1-preemptive tour servicing D, of length at most 0(log 2 n)\Q\ ■ B 
such that J2ieD Ti < 0(log n) • \D\B, where Tj denotes the total time spent in the vehicle by demand 
% G D. The bound on the delay uses the fact that max™ 1 d{si, U) < B. 

Choosing a large enough constant corresponding to p = ©(log n log m), the length of r is upper 
bounded by p\Q\ ■ B (since n < 2m). So the cutting procedure in Step ([3]) results in at most \Q\ pieces 
of r, each of length at most 2pB. The objects i G C" (as defined in Step Q) are called a cut objects. 
We restrict attention to the other objects C = D \ C" that are not 'cut'. For each object i G C', the 
path traced by it (under single vehicle tour r) from its source Sj to preemption-point and the path from 
its preemption-point to t{ are both completely contained in pieces of V. Figure [4] gives an example of 
objects in C and C" , and the cutting procedure. 

Claim 11 The expected number of objects in C" is at most Yli&D — 0( i O g m ) ' l-^l- 

rp 

Proof: The probability (over choice of r/) that object i G D is cut equals where Ti is the total time 
spent by i in tour r. The claim follows by linearity of expectation and J2i^D — 0(log n) ■ \D\B. ■ 

We can derandomize Step ([3]) and pick the best offset -q (there are at most polynomially many 



combinatorially distinct offsets). Claim 11 implies (again choosing large enough constant in p = 
0(lognlogm)) that \C'\ > (1 — 2I g TO )|D| > a ■ \D\ demands are not cut; note that this occurs 
with probability one since we use the best choice for m Now onwards we only consider the set C of 
uncut demands. Let V denote the pieces obtained by cutting r as above, recall \P\ < \Q\. A piece 
P G V is said to be non-trivial if the vehicle in the 1-preemptive tour r carries some C'-object while 
traversing P. Note that the number of non-trivial pieces in V is at most 2\C'\ < 2m: each C'-object 
appears in at most 2 pieces, one where it is moved from source to preemption-vertex and other from 
preemption- vertex to destination. Retain only the non-trivial pieces in V; so \V\ < min{|Q|, 2m} = z. 
The pieces in V may not be one-to-one assignable to the depots since the algorithm has not taken the 
depot locations into account. We determine which pieces may be assigned to depots by considering a 
matching problem between V and the depots in Step Q and ([6]). 



Load rebalancing. The bipartite graph H (defined in Step (|5]>) represents which pieces and depots 
may be assigned to each other. Piece P G V and depot / G Q are assignable iff d(f, P) < 2B, and in 
this case graph H contains an edge (P, /). We claim that corresponding to the 'maximal contracting' 
set S (defined in Step Q), the 2-matching ir (in Step Q) is guaranteed to exist. Note that |r(5)| < 
but |r(T)| > ^ for all T D S. For any T' C V \ S, let f(T') denote the neighborhood of T' in 
Q \ T(S). The maximality of S implies: for any non-empty T' <Z p\S, J|i + ^ = < 

\T(S UT')\ = \T(S)\ + |f (T')\, i.e. |f (T')\ > Hence by Hall's condition, there is a 2-matching 
7r : V \ S — > Q \ r(5). The set S and 2-matching it can be easily computed in polynomial time. 



Recursion. In Step ((7]), demands C' are further partitioned into two sets: C\ consists of objects that 
are either picked-up or dropped-off in some piece of S; and C2-objects are picked-up and dropped-off 
in pieces of V \ S. The vehicles T(5) suffice to serve all C\ objects, as shown below. 

Claim 12 There exists a schedule of vehicles T{S) serving demands C\, having makespan B. 

Proof: Consider the schedule of makespan B that serves all demands C = C\ U C2 using vehicles Q: 
this is implied by the promise on instance (Q, D, B). We claim that in this schedule, no vehicle from 
Q \ r(5) moves any C\ object. Suppose (for a contradiction) that the vehicle from depot / G Q \ T(S) 
moves object i G C\ at some point in this schedule; then it must be that d(f, Si) and d(f, U) < 2B. But 
since i G C\, at least one of Si or ti is in a piece of S, and this implies that there is some piece P G S 
with d(f, P) < 2B, i.e. / G T(5), which is a contradiction! Thus the only vehicles participating in the 
movement of C\ objects are T(S), which implies the claim. ■ 
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Figure 4: Cutting and patching steps in algorithm Partial. 



In the final schedule, a large fraction of C\ demands are served by vehicles r(<S), and all the C2 
demands are served by vehicles Q \ T(S). Figure [4] shows an example of this partition. 

Serving C\ demands. Based on Claim [l2| the recursive call Partial (T(S),Ci, B) (made in 
Step ([8])) satisfies the assumption required in Lemma 10 Since |r(«S)| < ^ < ^ < \Q\, we ob- 



tain inductively that Partial (r(<S), C\, B) returns a schedule of makespan (16 + 16p) • B covering at 
least a lsy • |Ci| demands of d, where y = min{|r(5)|, 2m}. Note that y < \T(S)\ < \V\/2 < z/2 
(recall \ V\ < z), which implies that at least a lgz_1 |Ci| demands are covered. 

Serving C2 demands. These are served by vehicles Q \ F(S) using the 2-matching ir, in two rounds 
as specified in Step ([8]). This suffices to serve all objects in C2 since for any i G C2, the paths traversed 
by object % under r, namely Sj ~» pi (its preemption-point) and pi -w ti are contained in pieces of V \ S. 
Furthermore, since |7r _1 (/)| < 2 for all / G Q \ T(S), the distance traveled by vehicle / in one round is 
at most 2 • 2(2B + 2pB). So the time taken by this schedule is at most 2 • 4(2B + 2pB) = (16 + 16p) • B. 

The schedule of vehicles T{S) (serving C\) and vehicles Q \ T(S) (serving Ci) can clearly be run 
in parallel. This takes time (16 + 16p) • B and covers in total at least |C 2 | +q lg2 ~ 1 |C l i| > a lgz_1 |C"| > 



a lgz |D| demands of D. This completes the proof of the inductive step of Lemma 10 



Using Lemma 10 repeatedly as mentioned earlier, we obtain an 0(log m ■ log n) -approximation 



algorithm. Using some preprocessing steps (described in the next subsection), we have Theorem[3] 



3.3 Weighted mDaR 

The multi- vehicle Dial-a-Ride problem as defined assumes that all objects have the same 'weight', i.e. 
each object occupies a unit capacity. In the weighted mDaR problem, each object i G [m] also has 
a weight Wi, and the capacity constraint requires that no vehicle carry a total weight of more than k 
at any time. In this section, we obtain an 0(log 3 n) -approximation algorithm for weighted mDaR. In 
particular this would imply Theorem [3] 

The algorithm for weighted mDaR first guesses the optimal makespan B of the given instance. It 
serves the demands in two phases: the 1st phase involves pre-processing 'heavy demands' and has a 



makespan of O(l) • B; the 2nd phase is identical to the algorithm Partial in Section 3.2 and covers all 
remaining demands. For every u, v G V let dem Ujt) denote the total weight of objects having source u 
and destination v. Define H = {(u, v) G V x V \ dem u v > k/2} to be the heavy vertex-pairs, and H 
the set of demands between pairs of H. 
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Phase I. In this pre-processing step, we cover all demands H between heavy vertex-pairs. The al- 
gorithm PreProc for solving instance U and its analysis follow along the lines of algorithm Partial in 
Section 3.2 in fact it is much easier. Below we use the same notation P C [q] for a set P of vehicles and 
the multi-set of depots corresponding to P. The algorithm PreProc will satisfy the following property. 



Lemma 13 Given any subset Q C [q] of vehicles and D C H of heavy vertex-pairs and bound B such 
that there is a schedule of vehicles Q covering demands D with makespan at most B, PreProc(0, D, B) 
returns a schedule of vehicles Q covering demands D that has makespan O(l) • B. 

Note: Invoking PreProc([g], H, B) gives the desired schedule covering H at makespan 0(1) • B. 

Algorithm PreProc first ensures that each edge of the minimum spanning tree on Q has length at 
most 3B, otherwise the problem decouples into two disjoint smaller problems (as in Step ([I]) of Partial). 
Hence, as in algorithm Partial, we also obtain that: 

The MST on all sources/destinations in D has length at most O(l) • \Q\B (1) 

Demands with source u and destination v are referred to as (u, v) demands. For every (u, v) G D, 
using a greedy procedure, one can partition all (u, v) demands such that the total weight of each part 
(except possibly the last) is between | and k. For any (u, v) G D, let g U}V denote the number of parts in 
this partition of (u, v) demands; note that g U)V < |dem Ujt) + 1 < Idem^ (since dem U) „ > k/2). The 
flow lower bound for the mDaR instance with vehicles Q and demands D equals Yl( u v)ed ^ em u,v • 
d(u, v) > jj^t J2( u v)ed 9u,v ■ d(u, v). Using the assumption in Lemma 
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we have that: 



9u,vd(u,v)<A\Q\B. (2) 

(u,v)eD 

In the rest of algorithm PreProc, we will consider each part in the above partition of (u, v) demands 
(for each (u, v) G D) as a single object of weight k; so all the objects in one part will be moved 
together. Scaling down the weights and capacity by k, we obtain the following equivalent unit-weight 
unit-capacity instance U: for each (u,v) G D there are g UjV demands with weight 1, source u and 
destination v, and each vehicle in Q has capacity 1. 

Since U has unit capacity and weights, we can use the Stacker-crane algorithm lfT6l to obtain a single 
vehicle non-preemptive tour t serving all demands, having length 1.8 times the Steiner and flow lower 
bounds. For this single vehicle problem, the lower bounds corresponding to U are: (Steiner bound) TSP 
on all sources/destinations in D, and (flow bound) J2(u v)eH 9u,v • d(u, v). From ([T| and ([2]) above, these 
are both 0(1) ■ \ Q\B. The algorithm next cuts tour r to obtain at most \Q\ pieces V such that each piece 
has length 0(1) • B. In addition it can be ensured that no object is being carried over the cut edges: 
since the vehicle carries at most one object at any time and each source-destination distance is at most 
B. When \Q\ = 1, the algorithm ends here, and returns this single tour of length 0(1) B that serves all 
demands D. This proves the base case of Lemma 
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When |0| > 2, construct bipartite graph % with vertex sets V and Q and an edge between piece 
P G V and depot / G Q iff d(f, P) < 2B; T(A) denotes the neighborhood of any A C V in H. The 
algorithm finds any maximal set S C V with |T(«S)| < |<S|/2 (as in Step ([5]) of Partial). This implies a 
2-matching ir : V \ S — ^ Q\ r(<S) such that there exists a schedule of vehicles T(S) serving demands 



in the pieces S with makespan B (c.f. Claim 12 1. Let CCD denote the heavy demand-pairs served 
in some piece of S. The final schedule returned by PreProc(0, D) involves: (i) schedule for vehicles 
Q \ r(<S) given by it (covering demand-pairs D \ C); and (ii) schedule for vehicles Y(S) obtained 
recursively PreProc(T(5), C) (covering demands C). The recursively obtained schedule has makespan 
O(l) • B by the induction hypothesis (since |r(<S)| < \Q\, and (F(S),C) satisfies the assumption in 



Lemma 13 1; hence the final schedule also has makespan 0(1) • B, which proves the inductive step. 
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Phase II. Let L = {(u, v ) G V x V | < dem u „ < A;/2} be the %/z? vertex-pairs. The algorithm 
for this phase treats all (u, v) demands as a single object of weight dem u ^ from u to v; s o th ere are 
m = \L\ < n 2 distinct objects. The algorithm for this case is identical to Partial of Section 3.2 for the 
unweighted case: Theorem [9] generalizes easily to the weighted case, and the Steiner and flow lower- 
bounds stay the same after combining demands in L. Thus we obtain an 0(log 2 m logn) = 0(log 3 n) 
approximate schedule that covers all remaining demands (setting m < n 2 ). 

Theorem 14 There is a randomized O (log 3 n)- approximation algorithm for weighted preemptive mDaR. 



3.4 Improved Guarantee for Metrics Excluding a Fixed Minor 

In this section, we show that our algorithm for preemptive mDaR achieves an 0(log 2 n) approximation 
guarantee when the metric is induced by a graph excluding any fixed minor. The main ingredient is the 
following improvement in the single vehicle structure from Theorem [9j 

Theorem 15 There is a randomized poly-time algorithm that for any instance of single vehicle dial-a- 
ride on metrics excluding a fixed minor, computes a 1 -preemptive tour r servicing all demands D that 
satifies the following conditions {where LB is the average of the Steiner and flow lower bounds): 

1. Total length: d{r) < O(logn) • LB. 

2. Bounded delay: J2ieD — 0(X) J2ieD <^( s i> U) where Tj is the total time spent by object i G D 
in the vehicle under the schedule given by r. 



Using this 1-preemptive tour within the algorithm of Section 3.2 (setting parameter p = 0(logm)), 



we immediately obtain Theorem |4] In proving Theorem [15} we first show that metrics induced by 



(fixed) minor- free graphs allow a so-called ^-separated cover (Subsection 3.4.1 1. Then we show (in 



Subsection 3.4.2 ) how this property can be used to obtain the single-vehicle tour claimed in Theorem 15 



3.4.1 7-Separated Covers 

We are given a metric (V, d) that is induced by a graph G = (V, E) and edge-lengths w : E — )■ Z + . 
For any pair u, v G V of vertices, the distance d(u, v) equals the shortest path between u and v under 
edge-lengths w. We assume that the graph G does not contain any i<C r -minor; here r is a fixed parameter. 
Recall that a cluster refers to any subset of vertices. We prove the following. 

Theorem 16 Given a K r -minor free graph G = (V, E, w) and integer 7 > 0, there is a polynomial- 
time algorithm to compute a collection C of clusters along with a partition Zi,...,Z p ofC such that: 

1. p = 2° {r \ 

2. For each I G [p] and distinct clusters A,B G Z\, the distance between A and B, d(A, B) > 7. 

3. The diameter of any cluster in S G C, max u ^ g g d(u, v) < 0(r 2 ) • 7. 

4. For any pair u,v 6 V of vertices with d(u, v) < 7, there is some cluster A G C having u,v G A. 
Such a collection C is called ^-separated cover. 

The proof of this theorem essentially follows from the KPR decomposition algorithm lF22l . In fact 
without conditions (1-2), it is implied by Theorem [7] Achieving conditions (1-2) requires a further 
modification to the KPR decomposition, as described below. 

We first remind the reader of the main property of the KPR decomposition [22J. We will consider 
graph G as having unit length edges, by subdividing each edge (u, v) G E to become a path of length 
w uv between u and v (this also increases the number of vertices). The distance function d remains 
unchanged by this modification. 
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Definition 17 (Theorem 4.2, |22|) LetG\ = GandS £ Z + a distance parameter. LetGi, G2, • ■ • , G r+ \ 
be any sequence of subgraphs of G, where each Gj+i is obtained from Gi as follows: 

1. Construct a breadth-first-search tree in Gifrom an arbitrary root-vertex. 

2. Select any set of 5 consecutive levels in this breadth-first-search tree, and let Gi+i be any con- 
nected component of the subgraph in G induced by these levels. 

Then, the diameter of 'G r +i is 0(r 2 ) • 5. 



We are now ready to proceed with the proof of Theorem 16 We give a recursive procedure to gen- 
erate the desired 7-separated cover C of G. We initialize C = 0, and invoke Split(G, V, 0, <fi). Algorithm 
Split takes as input: an induced subgraph H = (V(H), E{H)) of G, terminals T C V(H), depth of 
recursion i, and color r £ {0, 1, 2} 1 . Then it does the following. 

1. If i = r then add tuple (T, r) to collection C, and stop. 

2. Construct a BFS tree from any vertex of H, and enumerate the levels starting from zero (at the 
root). 

3. For each j £ {0, 1, 2}, integer I > 0, and connected component H' in the subgraph of G induced 
by levels {(3/ + j - 1)7, • ■ • , (31 + j + 3)7} of BFS do: 

(a) Set T «- vertices of T n V(H') in levels {(31 + j)7, • • • , (3/ + j + 2)7} of BFS. 

(b) Set t' <— t concatenated with j. 

(c) Recurse on Split(iJ', T' , i + 1, r'). 

At the end of the algorithm, C consists of several tuples of the form (C, r) where C is a cluster and 
r G {0, 1, 2} r is its color. This naturally corresponds to partitioning C into p = 3 r parts: clusters of the 



same color form a part. We now show that C satisfies the conditions in Theorem 16 Condition 1 clearly 
holds by the construction. 

Condition 4. Let u,v E V be such that d(u, v) < 7, and let P C V denote the vertices on any 
shortest u — v path in G (note that \P\ < 7 + 1). We claim inductively that for each i G {0, • • • , r}, 
there is some call to Split at depth i such that all vertices of P appear as terminals. The base case 
% = is obvious since all vertices are terminals. Assuming the claim to be true for depth i, we prove 
it for i + 1. Let (H, T, i, r) denote the call to Split at depth i where P C T C V(H). Since the P- 
vertices define a path of length at most 7 (in G and also H), all P-vertices appear in some 7 consecutive 
levels of any BFS on H. Thus there exist values j G {0, 1, 2} and I € Z + such that P is contained in 
levels {(3/ + ,7)7, • • • , (3/ + j + 2)7} of the BFS on H. This implies that one of the recursive calls 
(corresponding to this value of j and /) contains P as terminals. Using this claim when i = r implies 
that there is some cluster in C containing both u and v (in fact containing the entire path P). 

Condition 3. This is a direct consequence of the KPR decomposition. Note that every subgraph 



at depth r of algorithm Split is obtained from G via the sequence of operations in Definition 17 with 
S = 47 + 1. Thus each such subgraph has diameter 0(r 2 ) ■ 7. Finally, any cluster in C is a subset of 
some subgraph at depth r of Split; so it also has diameter 0(r 2 ) • 7. 

Condition 2. This well-separated property is the main reason for the use of 'terminals' in algorithm 
Split and the decomposition defined in Step ([3]). Let A, B be any two distinct clusters in C having the 
same color : we will show that d(A, B) > 7. Each of A and B corresponds to a "root-leaf path" in the 
recursion tree for algorithm Split (where the root is Split(G, V, 0, (ft)). Consider the call (H, T, i, r) to 
Split after which the paths for A and B diverge; let a a = (H a , T a , i + 1, r') and a& = (-H&, T&, i + 1, r') 
denote the respective recursive calls from (H, T,i,r) where A C T a and B C Tj (note that both have 
the same color r' since A and 73 have the same color at the end of the algorithm). We will show that 
d(T a , I},) > 7 which in particular implies d(A, B) > 7. Note that both a a and are generated by the 
same value j £ {0, 1, 2}, since they have the same color. Let a a (resp. a>b) correspond to value I = l a 
(resp. I = It). Consider the following cases: 
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1. l a 7^ If,. In this case, we show that the terminal-sets T a and are far apart in the BFS on 
H. Observe that T a appears within levels {(3/ a + • • • , (3/ a + j + 2)7}; and Tb within 
{(3Z&+ j)7, • • • , (3/;,+j + 2)7}. If Z a < (the other case is identical), then 3/fc+j > 3/ a + j + 3; 



i.e. T a and Tb are separated by at least 7 — 1 levels in the BFS of H. Using Claim 18 below, we 
have d(T a , Tb) > 7 in G; otherwise H contains a path from T a to Tb of length at most 7 — 1, 
meaning that T a and Tb are separated by < 7 — 2 levels in the BFS, a contradiction! 

2. l a = h = I- I n this case, it must be that H a and Hb are two disconnected components of the 
subgraph (of G) induced on levels {(31+ j— 1)7, • • • , (3/+j+3)7} of the BFS on H. Furthermore, 
T a and Tb both appear within levels {(3/ +,7)7, • • • , (3Z + j + 2)7}. Again by Claim 18 we must 
have d(T a , Tb) > 7 in G; otherwise H contains a path from T a to Tb of length at most 7 — 1, 
which would mean that H a and Hb are connected within levels {(31+ j — 1)7, • • • , (3/ + j + 3)7}. 

Claim 18 For every call Sp\'\t{H, T, i, r), and terminals x,y G T, ifd(x, y) < 7 (in graph G) then H 
contains an x — y path of length at most d(x, y). 

Proof: We proceed by induction on i. The claim is obviously true for the call Split(G, V, 0, cj>) at depth 
0. Assuming the claim for any depth i call Split(i7, T, i, r), we prove it for any call Split(i7', T' , i+1, t') 
generated by it. Let x, y G T' with d(x, y) < 7. Clearly x, y G T, and by the induction hypothesis, i? 
contains an x — y path ir of length at most d(x, y) < 7. It suffices to show that tt is also contained in H'. 
Let j G {0, 1, 2} and / G Z + denote the values that generated the recursive call Split(-£T, T', i + r'). 
Since x, y G T', both these vertices lie within levels {(3^ + .7)7, • • • ,(31 + j + 2)7} of the BFS on H. 
Furthermore, ir is an x — y path in H of length at most 7; so it is contained within levels {(3/ + j — 
1)7, • • • ,(3l+j + 3)7} and so it lies in the connected component H' that contains x and y. m 

Running time. It is easy to see that each call to Split generates only polynomial number of recursive 
calls. Since the depth of the recursion is at most r, the running time is polynomial for fixed r. 



This completes the proof of Theorem 16 



We note that the collection C also satisfies the 'sparse cover' property (i.e. condition (3) in Theo- 
rem|7]), namely each vertex v G V appears in at most 0(2 r ) clusters. It is easy to show (inductively) 
that the number of depth i calls to Split where any vertex v appears as a terminal is at most 2*. However 



this property is not required in the following proof of Theorem 15 



3.4.2 Algorithm for Theorem 15 



Recall that we are given an instance of single- vehicle Dial-a-Ride on a metric (V, d) induced by an 
edge-weighted AV-minor free graph G = (V, E, w); here w : E — > Z + denotes the edge-lengths and 
d the shortest-path distances under w. The vehicle has capacity k and {si, are the demand-pairs. 

Again r is a fixed constant. LB denotes the average of the Steiner (i.e. minimum TSP tour on all sources 
and destinations) and flow lower-bounds (i.e. I J2i^D d(si, U)) for this instance. 

Let A denote the diameter of the metric; by standard scaling arguments, we may assume WLOG 
that A is at most polynomial in n = \V\. We group the demands D into into [log 2 A] groups based 
on the source-destination distances: For each j = 1, • • • , [log 2 A] , group Gj consists of those demands 
i G D with 2 J_1 < d(si,ti) < 2 J . We will show how to service each group Gj separately at a cost of 
O(l) • LB such that the time spent by each object i G Gj in the vehicle is Ti = 0(2 J ) = 0(1) • d(s,i, ti). 
Since the number of groups is 0(log n), this would imply the theorem. 

Serving group Gj Set distance parameter 7 = 2 3 , and obtain a 7-separated cover C along with its 



partition {Z\, ■ ■ ■ , Z p }, using Theorem 16 Each demand i G Gj is assigned to some cluster of C 



containing both its end-points: since d(si,ti) < 7, this is well-defined by property Q of Theorem 15 



We further partition demands in Gj into Hi , • • • , H p where each Hi contains all demands assigned to 
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clusters of Z\. The algorithm will serve demands in each H\ separately by a tour of length 0(1) • LB 
that also satisfies the bounded delay condition 2. This suffices to prove the theorem since p = 0(1) by 
property ([T]) of Theorem 16 



Serving Hi. For any cluster S G Z\ let B(S) C Gj denote the demands assigned to S. Let LB(5) 
denote the average of the Steiner and flow lower bounds restricted to demands B(S). Any cluster S G Z\ 
is called non-trivial if B(S) ^ 0. Also pick an arbitrary vertex in each cluster as its center. Let T denote 
a 1.5-approximate TSP tour containing the centers of all non-trivial clusters of Zf, this can be computed 
in polynomial time ifTTTl . 

Claim 19 We have J2sez l LB ( S ) ^ ' LB > and d(T) < 0(1) • LB. 

Proof: It is clear that the sum of the flow lower-bounds over all clusters S G Z\ is at most the flow 
lower-bound on demands Hi C D. Let Tsp(5) denote the minimum length TSP on the end points of 
demands B(S), i.e. the Steiner lower bound on S. We show below that J2seZt "'"sp(S) < 0(1) ■ Tsp 
where Tsp denotes the minimum TSP on end points of Hi C D (i.e. Steiner lower bound on Hi). This 
would prove the first part of the claim. 

For any cluster S G Zi, let n(S) denote the number of connected segments in Tsp n S (i.e. Tsp 
restricted to S), and Of, • • • , ^(S) denote these segments. Note that tour Tsp travels from one cluster 



in Zi to another at least YlseZi n (S) ti mes - By property ([2]) of Theorem 15 the distance between any 
two distinct clusters in Zi is at least 7. So we obtain: 

Tsp > 7 n(S) (3) 
sez t 

We define a TSP tour on end-points of demands B(S) by connecting the end points of segments 
Cf , • • • , C^ s y This requires n(S) new edges, each of length at most 0(7) (i.e. diameter of clus- 



ter S, by property Q of Theorem 15 ). So the minimum TSP tour for any S G Zi has length Tsp(5) < 
E£? d(C s a )+ 0( 7 ) -n(S). Hence, 

n{S) 

Above, the second inequality follows from the fact that the Cj s are disjoint segments in Tsp, and 
the last inequality uses ([3]). 

To obtain the second part of the claim, augment Tsp as follows. For each non-trivial cluster S G Zi, 
add two edges to the center of S from some vertex in Tsp n S. This gives a TSP tour visiting the centers 
of all non-trivial clusters in Z\. Note that the number of non-trivial clusters of Zi is at most ^2 S( z Z[ n(S) 
since n(S) > 1 for all non-trivial S G Z\. Since the diameter of each S G Z\ is 0(7) (property Q of 
Theorem 15 I, the increase in length of Tsp is at most 2 • 0(7) X^SeZ; n (*^) — 0(1) ■ Tsp by Q. So 
there is a TSP tour on the centers of Zi's non-trivial clusters having length at most 0(1) • LB. ■ 

For each S G Zu the demands B(S) assigned to it are served separately. The flow lower-bound on 
S is at least ^\B(S)\ (recall that each demand i G Gj has d(si, U) > 7/2). Consider an instance T src 
of CVRP with all sources from B(S) and S"s center as the common destination. Since the diameter 
of S is 0(7), the flow lower bound of l src is at most 0(7) • < 0(1) • LB(5). The Steiner 

lower bound of I src is also O(l) • LB(S). So Theorem [8] (with delay parameter /3 = 2) implies a non- 
preemptive tour Tsrc(S) that moves all objects from sources of B(S) to the center of S, having length 
O(l) • LB(S') such that each object spends at most 0(1) • 7 time in the vehicle. Similarly, we can obtain 
non-preemptive tour Td es t(S) that moves all objects from the center of S to destinations of B(S), having 
length 0(1) • LB (S) such that each object spends at most 0(1) • 7 time in the vehicle. Concatenating 
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Tsrc(S) and Td est (S) gives a 1-preemptive tour serving demands in S having length O(l) • LB(5) such 
that each object i G B(S) spends at most 0(7) < 0(1) ■ d(si, U) time in the vehicle. 

The final solution for demands H\ traverses the TSP tour T on centers of all non-trivial clusters of 
Zi, and serves the demands in each cluster (as above) when it is visited. The resulting solution has length 
at most d(T) + 0(1) J2sez t LB ( 5 ') < 0(1) ' LB > bv Claim [l9| Moreover, each object i G Hi spends at 
most 0(1) • d(si, U) time in the vehicle. 

Finally, concatenating the tours serving each Hi (for I = 1,- • ■ ,p), and then the tours serving each 



Gj (for j = 1, • • ■ , [log 2 A] ), we obtain the 1-preemptive tour claimed in Theorem 15 
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